The four height variables, boundary correlations, and dissipative defects in the 
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We analyze the two-dimensional Abelian sandpile model, and demonstrate that the four height 
variables have different field identifications in the bulk, and along closed boundaries, but become 
identical, up to rescaling, along open boundaries. We consider two-point boundary correlations in 
detail, and discuss a number of complications that arise in the mapping from sandpile correlations 
to spanning tree correlations; the structure of our results suggests a conjecture that could greatly 
simplify future calculations. We find a number of three-point functions along closed boundaries, 
and propose closed boundary field identifications for the height variables. We analyze the effects of 
dissipative defect sites, at which the number of grains is not conserved, and show that dissipative 
defects along closed boundaries, and in the bulk, have no effect on any weakly allowed cluster vari- 
ables, or on their correlations. Along open boundaries, we find a particularly simple field structure; 
we calculate all n-point correlations, for any combinations of height variables and dissipative defect 
sites, and find that all heights and defects are represented by the same field operator. 

PACS numbers: 05.65.+b,45.70.-n 
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I. INTRODUCTION 

The Abelian sandpile model (ASM), introduced by 
Bak, Tang, and Wiesenfeld, is the original prototype for 
self-organized criticality [1]. Systems with self-organized 
criticality are naturally driven to a critical point, and 
thus can potentially explain how power laws occur in na- 
ture without any fine-tuning of parameters. Since their 
introduction, sandpile models have been used to model 
an extraordinarily wide range of systems, from earth- 
quakes [2] to river networks [3, 4]; see [5, 6] for reviews. 

To be precise, we are considering the two-dimensional 
isotropic Abelian sandpile model. This is a very simple 
model; in fact, its simplicity is its strength, since other- 
wise it could not act as a model for such a diverse range 
of physical systems. The ASM is defined on a square 
lattice, where each site has a height variable (the num- 
ber of grains of sand at that site) that can range from 
1 to 4. At each time step, a grain of sand is added to 
a random site. Any site with more than four grains is 
unstable, and collapses, losing four grains, and sending 
one grain to each of its neighbors. Unstable sites are 
repeatedly collapsed until all sites are stable. Then, a 
new time step begins — a grain is added to a random site, 
and the process begins anew [1] . initially, probabilities of 
configurations will depend on the initial conditions, but 
after a long period of time, the ASM develops a well- 
defined probability distribution of states, independent of 
the initial conditions [7] . Typically, the number of grains 
is conserved in each toppling, except for sites along open 
boundaries, where grains are lost with each toppling (i.e. 
fall off the edge). There must be at least one dissipative 
site — i.e. at least one site where the number of grains de- 
creases upon toppling — or else the sandpile would even- 
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tually reach a state where topplings continued endlessly 
during a single time step. 

Despite its simplicity, certain basic properties of the 
ASM remain unknown. For example, despite intensive 
work, the power law governing the sizes of avalanches 
in the ASM is still unknown — see [8] for a review. And 
while the height one variable is well understood, the roles 
played by the higher height variables (two, three, and 
four) are not. For example, no bulk two-point correlation 
functions of higher height variables are known. 

It is known that the ASM is related to the set of span- 
ning trees that can be drawn on the sandpile lattice, and 
that this relationship can be used to perform exact cal- 
culations of ASM probabilities [7, 9]. A spanning tree is 
a set of arrows drawn on the lattice, such that each site 
has exactly one arrow pointing from the site to a neigh- 
bor, and such that there are no closed loops of arrows. 
Following the path of arrows from any site will eventually 
lead off the edge of the sandpile (or, more generally, to 
a dissipative site, such as found on an open boundary) — 
the "site" off the edge of the sandpile is called the root. 
A number of relationships between the ASM and span- 
ning tree states are known. For example, the number 
of recurrent states of the ASM (states that occur with 
nonzero probability after a long amount of time) is equal 
to the number of spanning trees that can be drawn on 
the sandpile lattice [9]. 

Spanning trees are, in turn, related to the c = —2 
conformal field theory (CFT). The c = -2 CFT is the 
simplest known example of a logarithmic conformal field 
theory (LCFT), and is well understood [11-13]. 

A method introduced by Majumdar and Dhar ex- 
ploits the mapping between ASM states and spanning 
tree states to obtain exact ASM probabilities [14]. It 
has long been known that the Majumdar-Dhar method 
can be used to find the two-point correlation function 
of the unit height variable, which decays as 1/r 4 [14]. 
More recently, Mahieu and Ruelle used the Majumdar- 



2 



Dhar method to calculate correlation functions of a num- 
ber of ASM height configurations, known as weakly al- 
lowed cluster variables [15]. They not only found that 
all the correlations decayed as 1/r 4 , but were able to 
use their correlations to identify the thirteen simplest 
weakly allowed cluster variables with operators in the 
LCFT. These variables were all identified with linear 
combinations of three LCFT field operators, all of which 
had scaling dimension two, but only one of which — the 
8888 + 8888 operator — was isotropic. In some ways, this 
suggested that the higher height variables should be iden- 
tified with 8888 + 8888; on the other hand, Mahicu and 
Ruelle pointed out that this appeared inconsistent with 
LCFT operator product expansions (OPE's). 

Despite the power of these mappings, and of the 
Majumdar-Dhar method, fundamental questions about 
the ASM remain unanswered, because aspects of the 
mapping between the ASM and c = — 2 LCFT are still 
unknown — for example, it is not known what field oper- 
ators in the c = —2 LCFT represent the higher height 
variables of the ASM (or, indeed, whether such a rep- 
resentation even exists). A single site with height two, 
or any higher height variable, is not a weakly allowed 
cluster, and thus higher height probabilities and corre- 
lations cannot be calculated with the Majumdar-Dhar 
method. Priczzhev was able to extend the Majumdar- 
Dhar method to calculate the bulk probabilites for all 
higher height variables [10]. However, the bulk corre- 
lations of the higher height variables, which would be 
needed to obtain the field identifications of the higher 
height variables, remain unknown. 

Ivashkevich calculated all two-point correlation func- 
tions of all height variables, along open and closed bound- 
aries [16]. He found that all boundary correlations, be- 
tween all height variables, decayed as 1/r 4 , and argued 
that this implied that all four height variables should be 
represented by the same field operator (up to rescaling). 
Dhar has argued that, based on clustering properties of 
correlation functions, the bulk correlations should be ex- 
pected to factorize in a manner consistent with giving all 
four height variables the same field identification [17]. 

However, we argue here that the four height variables 
should in fact receive different field identifications, both 
along closed boundaries, and in the bulk, and propose 
field identifications along closed boundaries. Our con- 
clusions are based on analysis of closed boundary three- 
point functions, and of dissipative defect sites, as well as 
a reanalysis of the methods and results of Ivashkevich. 
However, we show that along open boundaries all four 
height variables, as well as dissipative defect sites, are 
represented by the same operator, 8988, in the c = — 2 
LCFT. We demonstrate this by computing all n-point 
correlations of height variables and dissipative defects. 

In sections II and III we briefly review the methods 
used by Majumdar, Dhar, and Priczzhev for studying the 
ASM. In section IV, we review Ivashkevich's calculations 
of the boundary height probabilities. 

In section V, and appendices A and B, we discuss issues 



associated with boundary correlation functions. While 
Ivashkevich has already calculated the boundary two- 
point correlations [16], we show that he mischaracterized 
the mapping between ASM configurations and spanning 
tree configurations, and a correct characterization results 
in a number of complications, necessitating a reanalysis 
of the two-point correlations. The relationship between 
ASM states and spanning tree states is not what one 
might have initially expected; we also note that linear 
relationships between nonlocal spanning tree conditions 
and local spanning tree conditions for one-point probabil- 
ities do not carry over in a simple fashion for multipoint 
correlation functions. Both of these complications intro- 
duce what we call "anomalous graphs" — while these com- 
plications are important, because they are technical in 
nature, we delegate much of the discussion to the appen- 
dices. In section V we calculate the anomalous graphs, 
and conjecture that the anomalous graphs have no ef- 
fect on the universal parts of any boundary correlation 
functions; while we have not been able to prove this con- 
jecture, it holds true for all correlation functions that we 
have calculated. 

In section VI, we look at correlation functions along 
closed boundaries. For two-point correlation functions, 
we find that while we disagree with Ivashkevich's rela- 
tionship between ASM and spanning tree states, we agree 
with his final results. However, we argue that these fi- 
nal results arc, in fact, not consistent with identifying all 
height variables with the same field operator. Next, we 
calculate all three-point functions along closed bound- 
aries that involve at least one unit height variable, and 
use these to make field identifications along closed bound- 
aries. Selected three-point functions appear in Eqs. (17- 
19), and we state the field identifications in Eqs. (20-22). 

Next, in section VII, we introduce the concept of a dis- 
sipative defect site, and discuss its effect on the lattice 
Green functions for the open, closed, and bulk cases. In 
section VIII, we show that in the closed and bulk cases, 
dissipative defects have no effects on any weakly allowed 
cluster variables. This demonstrates that an analysis of 
weakly allowed cluster variables, such as that in [15], can- 
not provide a complete picture of the ASM. Our results 
imply, as a particular case, that dissipative defects in the 
closed and bulk cases have no effect on the unit height 
probability, or on correlations of unit heights. They do, 
however, have an effect on the higher height variables; we 
show this analytically for the closed case, in Eqs. (28-29), 
and have checked this numerically for the bulk case. 

In sections IX, X, and XI, we compute all n-point cor- 
relation functions, for any number of height variables, 
and with any number of dissipative defects, along open 
boundaries. We find that there, all four height variables, 
and dissipative defects, are all represented by the same 
dimension two field, 8888. In fact, all local arrow dia- 
grams along open boundaries are represented by 8888, 
up to multiplicative prefactors. 

A short summary of these results can be found at [18]. 
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II. METHODS FOR ANALYSING THE ASM 

At its core, the ASM is a tractable model because the 
sandpilc model has an Abelian structure; the state of the 
sandpile does not depend on the order in which grains are 
added to the sites [7]. As a result of this Abelian struc- 
ture, it can be shown that the states of the sandpile fall 
into two simple categories. Some of the A N states of the 
sandpile (where N is the number of sites) are transient, 
which means that they can occur early in the ASM's evo- 
lution, but occur with zero probability after an infinitely 
long time. The other states are recurrent, and all occur 
with equal probability after long times. So the probabil- 
ity for a property X to occur is nothing more than the 
fraction of recurrent states having property X. 

To analyze the sandpile, it is convenient to allow more 
general toppling rules. We characterize the sandpile by 
a toppling matrix, A^j, where i and j are any lattice 

sites, i topples if its height is ever greater than A^, at 
which point its height goes down by Aj and the height 
of every other site j goes up by — A^j > (Ajj = if i 
and j are not neighbors). The original ASM, described 
in the introduction, has A^j = 4 when i = j (or Ajj = 3 
when i = j is along a closed boundary), Ajj = — 1 when 
i and j are nearest neighbors, and A^j = otherwise. 

Dhar was able to show that the number of recurrent 
states, given very general restrictions on A, is equal to 
dct(A) [7]. However, det(A) is also known to be equal 
to the number of spanning trees that can be drawn on 
the lattice [9]. In the spanning tree representation, Aj^ 
indicates the number of neighbors that the arrow from i 
can point to, A^j= —1 if an arrow can point from i to 

j, and A^ = otherwise. 

Certain height probabilities in the ASM can be equated 
with probabilities for spanning trees to have particular 
arrow configurations. Probabilities for some arrow con- 
figurations can be computed simply by modfiying the 
toppling matrix from A to A', in a way that enforces 
that arrow configuration. Then, the number of spanning 
trees with the configuration is det(A'), and the proba- 
bility of the configuration is det(A')/ dot (A). Defining 
B = A' — A, the probability becomes 



of sites on each end is decreased by one; so the three 
sites adjacent to i get maximum heights of 3, and i gets 
a maximum height of 4-3=1. It is not difficult to show 
that recurrent states S (of the original ASM) where i has 
height one are in one-to-one correspondence with the re- 
current states S' of the cut ASM. In this correspondence, 
we map from S to S' by lowering the heights of each of 
the three sites cut off from i by one. Letting j\, j 2 , and 
j3, be the three the neighbors that i has been cut off 
from, we have 
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Then the unit height probability is dct(I + BG) = 2(n — 
2)/7r 3 . This method was also used by Majumdar and 
Dhar to calculate the two-point correlation of the unit 
height variable [14]. 

Priczzhcv extended the Majumdar-Dhar method to al- 
low for the calculation of diagrams with closed loops. 
With the basic Majumdar-Dhar method, all off-diagonal 
entries of the toppling matrix are either or -1. Priez- 
zhev proved that if in A' we set n off-diagonal entries of 
A to — e, then 



lim 



det(A') 



(3) 



is equal to the number of arrow configurations such that 
each of the n corresponding arrows is in a closed loop 
of arrows, where each closed loop contributes a factor of 
— 1, and there are no closed loops other than those going 
through these n bonds. 

Such configurations are not spanning trees; spanning 
trees cannot have any closed loops. However, Priezzhev 
found that to calculate certain spanning tree probabili- 
ties, he needed to calculate graphs that had closed loops 
(^-graphs). We find this method useful for the calcula- 
tion of certain closed boundary correlations. 



dct(A') 
det(A) 



= det(I + BG). 



(1) 



G = A -1 is the well-known lattice Green function [24] 
(see appendix C). If A' only differs from A in a finite 
number of entries, then B is finite-dimensional, and the 
probability can be easily computed. 

Majumdar and Dhar used this method to find the prob- 
ability for a site i to have unit height [14]. To do this, 
they defined a modified, or "cut" ASM, in which three 
of the four bonds connecting i to nearest neighbors are 
removed. When a bond is removed, the maximum height 



III. HEIGHT PROBABILITIES 

Priezzhev determined a relationship between higher 
height probabilities and spanning tree states, which we 
review here [10]. 

Central to our analysis is the concept of forbidden sub- 
configurations (FSCs). A forbidden subconfiguration is a 
subset F of the lattice, such that for all i £ F, < c*(F), 
where is the height of site i, and c*(F) is number of 
neighbors that i has in F. Majumdar and Dhar proved 
that a state of the ASM is recurrent if and only if it has 
no FSCs [7, 9]. 
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The probability for a site i to have height two is more 
complicated than the height one probability [10]. In this 
case, changing the site height to one could either leave 
the ASM in an allowed (recurrent) state, or produce an 
FSC. The first case just gives the height one probability, 
which has already been calculated, so we consider the sec- 
ond case. Changing the height of i from two to one can 
produce multiple FSC's. Let F be the maximal forbid- 
den subconfiguration (MFSC) produced by this change. 
(Because more than one FSC can be produced, the word 
"maximal" is necessary for complete precision, and for 
this mapping to work; Priezzhev simply referred to "the" 
FSC, but this does not introduce any errors in his analy- 
sis [10].) F must contain i, and exactly one of the neigh- 
bors of i, and be simply connected, but can otherwise 
have arbitrary shape. The states S of the original ASM 
where changing the height of i from two to one produces 
F as the MFSC are in one-to-one correspondence with 
states S' of a modified ASM. In the modified ASM, all 
the bonds bordering F are removed, except for one (arbi- 
trarily chosen) bond of i. In the correspondence, we map 
from S to 5" by lowering heights of all sites that border 
F by the number of neighbors of F that they have been 
cut off from. (In this mapping, heights in F are unaf- 
fected.) With this mapping, the state S has no FSC's 
in the original ASM if and only if the state S' has no 
FSC's in the cut ASM. (Priezzhev's explanation used a 
slightly different, but equivalent, argument, based on the 
burning algorithm, a method for determining if a state is 
recurrent [7, 10]). 

The site % is called a predecessor of the site j in the 
spanning tree if the path from i to the root goes through 
j. We define NNPt-, as the number of nearest-neighbors 
of i that are predecessors of i. Then, the correspondence 
above shows that the number of states of the modified 
ASM is equal to the number of spanning trees of the 
modified lattice, which is in turn equal to the number of 
spanning trees where F is the set of predecessors of i. 
Summing over all possible sets F, we simply obtain the 
number of spanning trees where NNPt- = 1 . 

Similarly, it can be shown that the number of ASM 
states allowed when i has height h (or greater) but for- 
bidden when i has height h — 1 (or less) , is equal to the 
number of spanning trees where NNPt- = h— 1. Thus, the 
probability PasmQi) for the site to have exactly height 
h in the ASM is 
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where Ps P Tr(u — 1) is the probability that a randomly 
chosen spanning tree will have NNPt- = u — 1 , and rm is 

the maximum possible height of i. (m-* = 4 in the bulk, 
and along open boundaries, while toj = 3 along closed 
boundaries.) For more details, see [10]. 

This gives an exact representation of ASM height prob- 
abilities in terms of spanning tree probabilities. However, 





1 — i 


> — 1 







FIG. 1: Nonlocal arrow diagrams along closed boundaries. 



these spanning tree probabilities are not easy to calcu- 
late. Spanning tree probabilities that correspond to lo- 
cal restrictions on the spanning tree can be calculated 
with the Majumdar-Dhar method. However, the state- 
ment that NNPt- = u — 1 is a nonlocal restriction on the 
spanning tree (for u > 1). Priezzhev was able to calcu- 
late these nonlocal probabilities, but his calculations were 
complicated, and do not appear to be easily extensible to 
calculation of bulk correlations. However, this problem 
turns out to be more tractable along a boundary. 



IV. BOUNDARY HEIGHT PROBABILITIES 

For sites at the boundary, the relationship between 
height probabilities and NNP's still holds, and the NNP 
condition is still nonlocal. Nevertheless, Ivashkevich was 
able to show, through an ingenious transformation, that 
the ASM height probabilities are much easier to calculate 
along boundaries [16]. 

In figure 1 we list all possible nonlocal arrow config- 
urations around a site i of a closed boundary. In each 
picture, the dashed line is the boundary, and the central 
site is i. Large, solid, circles are predecessors of i, while 
large, open, circles are not. We see explicitly that the 
predecessor relationships are nonlocal. /3i and <j)\ differ 
only in whether the site above i leads to i by a chain of 
arrows — since the chain of arrows can go through sites 
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FIG. 2: Local arrow diagrams along closed boundaries. 



distant from i, this is a nonlocal distinction. If we can 
figure out the probabilities of all these diagrams, we can 
figure out the NNP probabilities (and thus the height 
probabilities). For example, the probability for NNPj- = 1 
is simply 2fa + 2fa + 27, since these diagrams catalogue 
all the ways that i can have exactly one NNP. 

These nonlocal diagrams are difficult to calculate. On 
the other hand, local restrictions are easily calculated 
with the Majumdar-Dhar method. All local arrow di- 
agrams along closed boundaries are shown in figure 2. 
Note that these diagrams do not have solid or open cir- 
cles, because predecessor relationships are not specified 
in local diagrams. Ivashkevich pointed out that the local 
arrow diagrams could be written as linear combinations 
of nonlocal arrow diagrams. For example, looking at fig- 
ures 1 and 2, we see that L Ci i = fa + fa. At first sight, 
there are more nonlocal arrow diagrams than local arrow 
diagrams, so such linear relationships would not appear 
to let us solve for the nonlocal arrow diagrams. However, 
Ivashkevich also pointed out that certain nonlocal arrow 
diagrams are equal in probability — for example, fa and 
fa are equal in probability, because we can make a one- 
to-one mapping from fa to fa by reversing all arrows 
in the long path of fa, and then switching the incom- 
ing arrow to i. Similarly, fa = fa. Then, we have as 
many nonlocal diagrams as local diagrams, and can solve 
for the nonlocal arrow diagrams. (In fact, along open 
boundaries, the number of local diagrams is one greater 
than the number of nonlocal arrow diagrams, so that the 
system is overconstraincd, providing a check on the cal- 
culations.) Ivashkevich used this to calculate all height 
probabilities along open and closed boundaries. See [16] 
for the full list of linear relationships between local and 
nonlocal diagrams. 
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FIG. 3: Anomalous graph of the first kind arising in the cal- 
culation of the two-point function. 



V. BOUNDARY TWO-POINT CORRELATIONS 
AND ANOMALOUS GRAPHS 

The calculation of boundary correlations is much more 
difficult. We show in appendix A that Ivashkevich's cal- 
culation of the two-point functions was incorrect, and 
ignored complications that arise in the relationship be- 
tween ASM height correlations and spanning tree correla- 
tions (although his end result turns out to be correct). In 
appendix B, we discuss further complications that arise 
in transforming from nonlocal spanning tree correlations 
to local spanning tree correlations. We summarize the re- 
sults here, and analyze the resulting "anomalous graphs" . 

The first complication arises in the correspondence be- 
tween ASM height probabilities and spanning tree prob- 
abilities. It would be natural to think that, analogously 
to Eq. (4), the ASM probability, PASM(h^,hj), for the 

sites i and j to have exactly heights and hj should be 
given by 
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(m| + 1 — u) (mj + 1 — v) 



(5) 

where Ps P Tr(u — ljf — 1) is the probability that in a 
spanning tree, NNPj- = u— 1 and NNPj- = v — 1. However, 

this turns out to not be quite the case. Eq. (5) is a natural 
guess, which we call a "naive" approach, but as shown in 
appendix A, the left and right sides of Eq. (5) differ by a 
subset of spanning trees that we call anomalous graphs 
of the first kind. (These graphs are not anomalous in 
any physical sense; we simply mean that they differ from 
what we would get, using a certain naive starting point.) 

The set of anomalous graphs of the first kind, for the 
closed case, is shown in figure 3. (We represent the root 
with a star.) In the graph, NNP r = NNPj = 1, so this 
graph appears in the right-hand side of Eq. (5) for = 
hj = 2. However, we show in appendix A that this graph 
does not contribute to the 2-2 (height two-height two) 
correlation, but instead contributes to the 2-3 and 3-2 
correlations, and gets subtracted from the 3-3 correlation. 

Second, leaving aside for now the anomalous graphs of 
the first kind, we need to calculate correlations of nonlo- 
cal arrow diagrams. It would be convenient if we could 
use the linear relationships relating nonlocal arrow dia- 
grams to local arrow diagrams found for one-point func- 
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FIG. 4: Anomalous graphs of the second kind arising in the 
calculation of the two-point functions. 
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FIG. 5: Anomalous graph & as a linear combination of closed 
loop diagrams. 
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x is the separation between i and j along the defect. 
g(x) is the Green function between i and j, and diverges 
as ln(L), where L is the system size (it also diverges 
as ln(x)). The restriction that spanning trees should 
have no closed loops greatly limits the number of pos- 



tions (section IV), and use them independently at i and 
at j for two-point functions. We again call this approach 
"naive," and again, this approach does not quite work. 
The problem arises because for one-point functions, we 
treated 4>\ and (f>2 as equivalent, based on a one-to-one 
correspondence in which a long path was reversed. In a 
correlation function of nonlocal arrow diagrams, the long 
path from a (f>i at i may go through arrow constraints 
near j, which are not free to be reversed. We discuss this 
problem in detail in appendix B. Consideration of this 
problem shows that, relative to the naive approach, our 
results are changed by graphs a and 6, shown in figure 4. 
We call these anomalous graphs of the second kind. 



The anomalous graphs a, 6, and c can be calculated 
with the extension of the method of Priczzhcv, discussed 
in section II [10]. We discuss only the calculation of the 
6 term; the analysis of the other terms is similar. 

6 represents a subset of spanning trees, and thus cannot 
have any closed loops. However, it comes "very close" to 
having a closed loop that includes the distant sites i and 
j, and we see in figure 5 that b can be written as a sum 
of closed loop diagrams. 



Priczzhev's method allows us to calculate the closed 
loop diagrams. We represent an arrow whose weight in A 
is set to — e (e — > oo) with a wavy bond line. As discussed 
in section II, these bonds must be part of a closed loop, 
and we get a factor of —1 for every closed loop. This 
gives the relations in figure 6. Taking the difference of 
the two graphs in figure 6 then gives the value of a closed 
loop diagrams that goes through both i and j. Using this 
method, we find the number of diagrams bl, 62, and 63 
(labeled in figure 5), as ratios of N, the total number of 
spanning trees: 



sible spanning trees, when the outlets to the root (open 
boundaries) are very far away. So diagrams such as 61, 
62, and 63, that allow a closed loop, are much more nu- 
merous than diagrams of spanning trees. 

However, to find 6, we take the linear combination, 



(3tt - 8) 2 (tt. 9 (x) - 1) (-128 + 48tt + tt 2 ) + (256 - 192tt + 307r 2 )7r 5 (a;) 1 

(3tt-8) 2 (^(x)-1) (3n-8)(2ng(x)-l) (3tt - 8)((4 - tt) + (3tt - 8)ng(x)) 1 

2tt 4 x 2 4tt 3 x 3 + 27r 4 a; 4 + { x 5> [) 

(3tt - 8) 2 (7r. 9 (a;) - 1) (3tt - 8){2ng(x) - 1) (-128 + 48tt - tt 2 ) + (256 - 192tt + 427r 2 )7r 5 (x) 1 

4^2 + + + y-tf) 

(8) 



_ SSL _ = 

= © rO 

FIG. 6: Use of -e weight bonds to evaluate closed loop diagrams. 



(Nbi — Nb2 + Nb3)/N, and the ln(L) divergences cancel 
(this provides a check on our calculations), a and c can 
be found similarly. We find 



VI. TWO- AND THREE-POINT CLOSED 
BOUNDARY CORRELATION FUNCTIONS 

We define, for all correlation functions along closed 
boundaries, 



a = 



b = 



c = 
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a and 6 are both of order 1/x 4 . The two-point correla- 
tion functions turn out to decay as l/x 4 , so the anoma- 
lous graphs could, in principle, affect the universal parts 
of the correlation functions. However, the anomalous 
graphs of the second kind come in the combination (a — b) 
(see Eqs. (B5-B7)). So their total contributions to the 
two-point correlations are 0(l/x 5 ), and can be dropped. 

The end result is rather surprising. A "naive" ap- 
proach might simply apply the relationship between ASM 
states and NNP conditions found for the one-point func- 
tions, independently at i and j (i.e. extend Eq. (4) to 
Eq. (5)), and then apply the relationships between non- 
local arrow diagrams and local arrow diagrams found for 
the one-point functions, independently at i and j. Nei- 
ther of these steps is correct, and a correct analysis pro- 
duces correction terms (the anomalous graphs) to this 
naive approach. But, somehow, the anomalous graphs, 
while nonzero, produce no correction to the leading- 
order, universal results at any stage of the computation; 
the naive approach gives the answers. In fact, we find in 
the following sections that the naive approach again gives 
correct results for all three-point closed boundary correla- 
tions that we have calculated, and for all open boundary 
correlations. This leads us to conjecture that the naive 
approach always produces correct universal results, for 
all correlations. If this conjecture were proven true, it 
would greatly simplify further calculations — for example, 
the anomalous diagrams have prevented us from calcu- 
lating the 2-2-2 correlation along closed boundaries. 



/ c (ai,a 2 



,a n ) = 

Poi.c) • 



Pa„,c))c , (12) 



In this correlation function, the height h Xu at the bound- 
ary site x u is required to be a u . We have subtracted off 
the constant boundary probabilities, p au .c, which were 
found in [16], as described in section IV. The subscript 
"c" stands for "closed." As already noted, despite errors 
in the setup in [16], the results of [16] are nevertheless 
correct, where it was found that 



/c(M) 

/c(l,2) 
/c(2,2) 



9 48 64 
12 68 96 



7T Z 



7T° 



7P 



1 



{xi -x 2 y 



(14) 



61 96 144\ 
+ ^ I 



4ir 2 7T 3 7r 4 / (xi — X2) 4 



+ ... 



(15) 



The correlation functions involving the height three vari- 
ables were also calculated, but we do not write them here, 
as they are determined by the requirement that all height 
probabilities must sum to one at every site. (There is a 
misprint in the result for / c (3,3) in [16].) 

Ivashkevich argued that the fact that all two-point cor- 
relations decay as l/(xi — X2) 4 indicates that all three 
height variables are represented by the same field opera- 
tor. However, if all height variables were represented by 
the same operator, we would expect the two-point func- 
tions to factorize, as 



fc(u,v) 



(xi - x 2 ) 4 



u,vG {1,2,3}, (16) 



for some constants K u . However, the results in Eqs. (13- 
15) do not factorize in this manner. Dhar argued that 
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we should expect this factorization for bulk correlations, 
based on the "clustering properties of correlation func- 
tions," but we see that this factorization already fails 
along closed boundaries [17]. (We will see later that the 
open boundary correlations do, however, factorize in this 
manner, for all n-point correlations.) 

To clarify the field identifications, we have calculated 
all three-point functions along closed boundaries, where 
at least one of the heights is the unit height. Some of our 
results are 



/c(l,l,l) 

/c(l,l,2) 
/c(l,2,2) 



2(3tt-8) s 



7r 6 (xi - x 2 ) 2 (x 1 - x 3 ) 2 (x 2 - x 3 ) 2 

8(7T-3)(37T-8) 2 
ir e (xi - X2) 2 {XI - X 3 ) 2 (x 2 - X3) 2 

(37T-8) 2 ^ 
7r 5 (a;i — x 3 ) 3 (x2 — X3) 3 

4(37T-8)(-57r 2 + 397r-72) 

TT & (Xl - X 2 ) 2 {XI - X 3 ) 2 (x 2 - X 3 ) 2 

(37r- 8)(24 - 77r) 



+ 



2ir 5 (xi - x 2 ) 3 (xi - x 3 ) 3 



+ ... 



(17) 



(18) 



(19) 



Other three-point correlation functions, calculated 
with the same methods, are listed in appendix D. They 
are all consistent with the requirement that the three 
height probabilities must sum to one at any site, and 
permutation symmetry, thus providing a check on our 
calculations. 

We again get a number of anomalous graphs (relative 
to a naive approach), and as stated in the previous sec- 
tion, again find that all anomalous graphs cancel in the 
universal, leading-order terms of the correlation function. 

These correlation functions are consistent with identi- 
fying the height variables with the following field opera- 
tors in the c = -2 CFT: 



Height one : 
Height two : 



Height three : 



2(3tt- 



7T 

6(tt -4) 



(8686) 



(8686) + —68 2 6 



(8686) 



2tt 



2tt 
-68 2 6 



(20) 
(21) 
(22) 



The representation of the c = — 2 CFT used here is de- 
scribed briefly in appendix E. Note that the boundary 
correlations in Eqs. (17-19) are the same as the bulk cor- 
relations of Eqs. (20-22), and that while the c = -2 CFT 
contains holomorphic and antiholomorphic fields (the 8 
and 5 of Eq. (El)), the fields in Eqs. (20-22) contain 
only holomorphic fields. This is consistent with boundary 
CFT. While fields in the bulk generally have holomorphic 
and antiholomorphic parts, near a boundary the antiholo- 
morphic pieces behave, in all correlation functions, like 
holomorphic pieces at mirror positions across the bound- 
ary [19]. 



It is also consistent to make the substitution 6^8, 
8 — > —6 in these field identifications, as the c = —2 LCFT 
is symmetric under this transformation (see Eq. (El)). 

The fact that the field identifications for the height 
variables differ along a closed boundary proves that they 
must also differ in the bulk. This is because in a CFT, 
boundary operators are derived from operator product 
expansions (OPE's) of bulk operators [19]. Furthermore, 
in appendix F we present a simple argument, based on 
general CFT principles, and not on any detailed calcula- 
tions, that the height variables must have different field 
identifications in the bulk. 

We have not been able to calculate three-point corre- 
lation functions that have no unit height variables. The 
basic problem is with the anomalous diagrams that arise 
when we convert from nonlocal arrow diagrams to local 
arrow diagrams (as in appendix B). The trick shown 
in figure 6, for evaluating the resultant closed loop dia- 
grams, does not work for these three-point functions. We 
note that if we use the conjecture proposed in section V 
(i.e., ignore the anomalous graphs), we obtain 



/ c (2,2,2) = - 



(24 - 5tt)(-576 + 384tt - 61tt 2 ) 
47r 6 (xi - x 2 ) 2 (x 1 - x 3 ) 2 (x 2 - x 3 ) : 



+ 



(23) 

(and other three-point functions consistent with the re- 
quirement that all three height probabilities must sum 
to one at any site). This correlation function is consis- 
tent with the field identification in Eq. (21), providing 
support for our conjecture. 



VII. DISSIPATIVE DEFECT SITES, GENERAL 

We now consider the effects of dissipative defects on 
the ASM. Generally, at sites in the bulk, or along closed 
boundaries, the number of grains is conserved at each 
toppling. Usually, it is only at open boundaries that 
the number of grains is not conserved; there, of the four 
grains toppled, three grains are sent to neighbors, while 
the remaining grain goes off the edge of the sandpile, to 
the root. 

Some dissipation (i.e. sites where topplings remove 
grains from the sandpile) is necessary for the sandpile 
model to be well defined, since otherwise we would end up 
with states where the topplings never terminated. Nev- 
ertheless, dissipation often plays a minor role in analysis 
of the sandpile, because properties are often studied in 
the bulk of the ASM, with the dissipative sites along the 
open boundaries infinitely far away. 

Some previous studies have investigated the effect of 
adding dissipation throughout the bulk of the ASM. In- 
stead of having bulk sites topple when their height is 
greater than 4, they topple when their height is greater 
than 4 + k (k > 0). Then, one grain is sent to each 
of the four neighbors, and k grains are lost to the root. 
It has been shown, both numerically and analytically, 
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that if this is done at all sites, the ASM is taken off 
the critical point, and the power law correlations are de- 
stroyed [20-22]. This happens even when k is infinites- 
imal. (Although this modification to the ASM has its 
most obvious interpretation for integer k, the theory can 
be given a sensible interpretation for any rational value 
of k. See [20] for details.) More recently, Mahicu and 
Ruelle have demonstrated the precise manner in which 
dissipation throughout the bulk takes the ASM off the 
critical point. They found that the dissipation has ex- 
actly the same effect on correlation functions of weakly 
allowed cluster variables, as adding the integral of the di- 
mension variable, 99, to the c = — 2 CFT [15]. Adding 
dissipation along a line has been shown to split the ASM 
into two separate half-planes, each with open boundary 
conditions [23]. 

Here, we consider the effect of adding dissipation at 
only a single defect site. Then, the methods of Majumdar 
and Dhar still work, but we need to use a modified lattice 
Green function. If k grains of sand are dissipated at the 
lattice position d, then we call k the "strength" of the 
defect. The toppling matrix is then changed from the 
defect- free toppling matrix, A , to 

The Green function is simply the inverse of the toppling 
matrix, and is changed from the defect-free Green func- 
tion, Go (described in appendix C), to 

G(i, j) = G (i, J) 1 * G Q {i, d)G (d, j) (25) 
1 + kG (d,d) 

This holds for any value of k, and regardless of the loca- 
tion of the defect. Nevertheless, the defect behaves very 
differently in the open case, and in the closed/bulk cases. 
This is because the Green function between nearby lat- 
tice sites is 0(1) near an open boundary, but O(lnL) 
near a closed boundary, or in the bulk [24, 25]. L is the 
system size, or more generally, is of the same order-of- 
magnitude as the distance to the nearest open boundary. 
This divergence in the Green function as L — > oo for the 
closed and bulk cases is usually not an issue, since in most 
cases, we are concerned with differences in Green func- 
tions. However, here the divergence of all the Go terms 
makes Eq. (25) unwieldy, although technically correct. 
(Eq. (25) can be used in the open case without modifica- 
tion.) Wc work in a limit where the distances between i, 
j, and d, while possibly large, are all much less than L. 
In this limit, dropping terms of order l/(lnL), Eq. (25) 
becomes 

G(i, j) = Go (I J) - Go (i, d) - Go (d, j) + Go (d, d) (26) 

Note that Eq. (26) is independent of k. This makes sense, 
since in the bulk, or along a closed boundary, spanning 



trees have to travel far to reach the root. But with the 
defect given by Eq. (25), k bonds are added from the 
defect d to the root. Adding a dissipative defect provides 
such an "easy" way to reach the root, that with high 
probability (probability one as L — > oo), all nearby points 
will be predecessors of the dissipative defect, regardless 
of the value of k. The set of spanning trees will thus 
be the same, in the L — > oo limit, for any k. Note also 
that the Green function in Eq. (26) no longer diverges as 
L — > oo, which is appropriate, as we are no longer O(L) 
from any dissipative sites. 

VIII. DISSIPATIVE DEFECT SITES, CLOSED 
AND BULK CASES 

Surprisingly, it turns out that a dissipative defect, ei- 
ther in the bulk, or on or near a closed boundary, has 
no effect on any weakly allowed cluster variables in the 
ASM. Weakly allowed cluster variables are height config- 
urations that result in a subconfiguration that contains 
an FSC if any height in the configuration is reduced by 
one [26]. Examples of weakly allowed cluster variables 
are a single height one variable, or a height one adjacent 
to a height two. Such variables can be calculated with 
the Majumdar-Dhar method by the removal of a set of 
bonds in the ASM/spanning tree. We note that corre- 
lations of weakly allowed cluster variables (such as all 
correlations of the unit height variable) are also weakly 
allowed cluster variables. 

Probabilities of weakly allowed cluster variables can be 
calculated as det(I + BG), as in section II. To analyze 
the effects of the defect, we want to consider the effect of 
modifying the Green function from the defect-free Green 
function Go, to the Green function in Eq. (26), for a 
fixed matrix B (i.e. for a specific weakly allowed cluster 
variable). 

In general, for local arrow restrictions, each row of B 
must sum to zero, because if the restrictions on the span- 
ning trees prevent an arrow from i to j, then 5jj goes 
down by 1, while B^j goes up by 1. (For example, for the 
height one variable, the matrix B in Eq. (2) arises from 
the restriction that no arrows can point from i to j\, 
or js, nor from j\, or j'3 to i.) For the weakly allowed 
cluster variables, B is symmetric, since if the arrow from 
i to j is forbidden, then so is the arrow from j to i. So 
every column of B also sums to zero. 

Since every row of B sums to zero, the parts of G that 
are independent of the row index of G make no contribu- 
tion to BG, and thus no contribution to the probability 
dct(I + BG). And det(I + BG) = dct(I + GB), so since 
every column of B sums to zero, the parts of G that are 
independent of the column index also make no contribu- 
tion to the probability. The last three terms of Eq. (26) 
all depend either only on the row index, or only on the 
column index. So a dissipative defect has no effect on any 
weakly allowed cluster probabilities (cither on one-point 
probabilities or on correlations). 
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As a special case, this means that the unit height prob- 
ability, and its correlations, are unaffected by closed or 
bulk dissipative defects. However, the higher height vari- 
ables are affected. Using the Green function in Eq. (26), 
and the methods described in section IV, we find that 
along a closed boundary, with a defect at the origin, we 
have the following height probabilites at xf. 



/c(l) = (27) 



/c(3) = +^ ! + ... (29) 

We have numerically confirmed these results. These re- 
sults provide further evidence that the height two and 
three variables have different field identifications along 
closed boundaries. 

Since the height two and three variables have dimen- 
sion two, this indicates that a dissipative defect along a 
closed boundary is a dimension zero operator. Consis- 
tent with this, uniform dissipation in the bulk has been 
identified with the integral of the dimension zero oper- 
ator 99 [15]. However, the correlation of 99 with the 
height two and three operators in Eqs. (21-22) does not 
produce the correlations in Eqs. (27-29); this situation 
requires further analysis. 

In the bulk, we would also expect that the higher 
height probabilites would be affected by a defect site, 
and have confirmed this with numerical simulations, al- 
though have not proven this analytically. 

The fact that weakly allowed cluster variables have no 
correlations with bulk or closed defects provides com- 
pelling evidence that weakly allowed cluster variables do 
not provide a complete picture of the sandpile model. 
This has particular bearing on the analysis of Mahicu 
and Ruellc [15]. They studied specific bulk correlations 
of the simplest weakly allowed cluster variables, and de- 
veloped a complete field picture for these variables. They 
found that (at the critical point) these variables are all 
linear combinations of three dimension two variables, 
8989 + 8989, 8989, and 8989, strongly indicating that all 
weakly allowed cluster variables are linear combinations 
of these three fields. However, this analysis left the status 
of the height two variable unresolved. Mahieu and Ruellc 
pointed out that since the height two variable appears in 
a number of the weakly allowed cluster variables, it might 
be expected that the height two variable would also be a 
linear combination of these three fields, or more specifi- 
cally, proportional to the sole rotationally invariant field, 
8989 + 8989 [15]. But they also noted that such an iden- 
tification appeared inconsistent with the fusion rules of 
the c = — 2 CFT, which would indicate a different field 
identification. The analysis here points strongly to the 
latter conclusion, although the specific field identification 
in the bulk remains unresolved. 



IX. ALL n-POINT CORRELATIONS ALONG 
OPEN BOUNDARIES, PART I 

We have calculated all n-point correlations of all four 
height variables, along open boundaries, in the presence 
of an arbitrary number of dissipative defects. We begin 
by discussing why this case is so tractable (in contrast to 
the closed case, where we have been unable to calculate 
the three-point function of the height two variable). 

The heights of the correlation function are placed at 
x\, X2, ■ ■ ■ , x n , and defining x a b = x a —xt, = c^x, we work 
in the limit x — > oo, where the c a t's are kept constant. 

As discussed in section V, and appendices A and B, 
a number of anomalous terms arise in the computation 
of correlation functions. While the discussion in these 
sections focused on closed boundary correlations, similar 
anomalous graphs arise in open boundary correlations. 
However, it turns out that these anomalous graphs pro- 
duce no contributions to the universal parts of any cor- 
relation functions, greatly simplifying matters. We prove 
this claim in this section, and in the next section look at 
the actual calculation of the correlation functions. 

We start by focusing on the two-point correlations. 
Note that the anomalous graphs found thus far, in fig- 
ure 3 and 4, all involve "nearly-closed" loops: the trees 
have paths that go from the neighborhood of i to the 
neighborhood of j, and from the neighborhood of j to 
the neighborhood of i. The paths do not actually form 
closed loops, since no closed loops are allowed in spanning 
trees, but they do come very close (within one site) . The 
reasons for this are general, so similar structures will arise 
in all anomalous graphs, for all correlation functions. For 
example, the anomalous graphs in figure 4 arose because 
a long, nearly-closed loop from one site could not be re- 
versed in direction, if it passed through fixed arrows at 
the other site (see figure 10 of appendix B). 

In the open case, these anomalous graphs between i 
and j always fall off faster than 0{l/x A ). This is in 
contrast to the closed boundary case, where such dia- 
grams diverge — see Eqs. (6-8). The difference results 
from the Green functions. While the Green function 
diverges as ln(x) along closed boundaries, it decays as 
1/x 2 along open boundaries (see appendix C). Using 
Priczzhev's method, the matrix determinant for evalu- 
ating any closed loop diagrams necessarily involves two 
Green functions, one from i to j, and another from j 
to i, giving an overall factor of 1/x 4 . Furthermore, cal- 
culating the diagrams requires two matrix determinants, 
which come with leading terms equal in magnitude, but 
opposite in sign — see figure 6. The 0{l/x A ) parts of the 
closed loop diagrams thus cancel along open boundaries. 
So the anomalous graphs for the two-point functions au- 
tomatically fall off faster than 0{l/x A ), and do not need 
to be considered when calculating leading-order, univer- 
sal parts of correlation functions. 

By this logic, for any n-point open boundary correla- 
tions, any anomalous graphs must have terms that decay 
as l/(x a — Xb) p , where p > 5, for some a, b G 1, 2, • • • n. 
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Aside from the sites at x a and Xb, there are (n — 2) other 
sites that need to appear in the connected correlation 
function. Each brings a new Green function, of 0(l/x 2 ), 
so the overall contribution of any anomalous graph must 
decay at least as fast as £>(l/a; 5+2 ( n - 2 )) = 0(l/x 2n+1 ). 
But we will see in the next section that all n-point cor- 
relations decay to leading order as 1 /x 2n . So the anoma- 
lous graphs have no effect on the universal parts of any 
n-point correlation functions. The conjecture at the end 
of section V has thus been proven for all open boundary 
correlations. 



X. ALL n-POINT CORRELATIONS ALONG 
OPEN BOUNDARIES, PART II 

Since we can ignore the anomalous graphs for open 
boundary correlation functions, no error is introduced 
by writing the height probability at each site as a linear 
combination of local arrow diagrams, independently us- 
ing at each site the linear relationships derived for the 
one-point functions. Defining the open boundary corre- 
lation / op analogously to f c for the closed case (Eq. (12)), 
we then have 



Moc N loc N loc 
— ^ ' ^ ' ■ • • ^ ' D aiUl D a2U2 . . . Da n u n 
ui—1 U2 — 1 u n —l 

(x n )) (30) 

Each L oPtUf (xf) represents a local arrow diagram at Xf, 
analogous to the diagrams in figure 2, but for the open 
case, and Ni oc is the total number of possible local ar- 
row diagrams at a single site (see [16] for the list of dia- 
grams). D is a constant matrix expressing height prob- 
abilities in terms of local arrow diagrams, for one-point 
functions, and was (implicitly) found in [16]. Each cor- 
relation of local arrow diagrams can now be calculated 
with the Majumdar-Dhar method. 

If a site i has local arrow constraints u, we express 
those constraints by a matrix B u , and let G uu be the 
Green function matrix for the sites around i. B„ and G uu 
are both associated only with sites in the vicinity of i. 
p u = dct(I + B„G utI ) gives the one-point probability for 
the local arrow diag ram Lo-p.w The two-point correlation 
of local arrow diagrams ui and Ui is given by det(I+BG), 
where B is block diagonal, with B Ul and B„ 2 along the 
block diagonal, and G is made of the four matrix blocks 



G Ul Ul , G^ 



Gu 2 ui , and G. U2 



Mahicu and Rucllc 



/o P (ai, «2, ■ • ■ , a n ) 



found that the leading order contribution to the bulk two- 
point probability is given by [15] 



dct(I + BG) = Trace • 



I 



-B M , C 



I 



-B M „ G„ 



1 + D U2 G, 



t H ~t~ B Ul G UlUl ii "T" *->u 2 ^*U 2 U 2 

Similarly, they found that the bulk, leading-order, contribution to the three-point probability is given by 



(31) 



dct(I + BG) = 

PuiPu 2 Pu 3 Trace 

p Ul Pu 2 Pu 3 Tia,ce 



I 



i 



I ~T~ Gu2U2 

I 



-B U2 G U2U3 : 



B U3 G« 3M3 



G U3Ul 



'I + B U3 G 



I 



! I + B U2 G 



U2 ^*U2U2 



(32) 



(Eqs. (31-32) are written in a different form than the 
expressions in [15], but are equivalent.) 

The derivation in [15] of Eq. (31) in the bulk relied 
on the fact that the leading-order contribution to the 
two-point function comes from the pieces of det(I + BG) 
with two terms off the block diagonal (i.e. one term from 
G UlU2 , and one term from G U2Ul ). Similarly, the deriva- 
tion of Eq. (32) was based on the fact that the leading- 
order, connected, contribution to the three-point func- 
tion comes from the terms of det(I + BG) with three 
terms off the block diagonal. 

The trace formulae can be extended for all higher- 



order correlations for the open case. We will see that the 
leading-order contribution to the open boundary n-point 
function decays as 0{l/x 2n ). The open boundary Green 
function (appendix C) between (x\,yi) and (£2,2/2) i s 

n i \ (yi + 1 )(y2 + i) . 

G op ,o(a;i,yi; 3:2,2/2) = -, ^ h... (33 

7r(xi - x 2 y 

Here, x labels distance along the boundary, and y labels 
distance from the boundary (the boundary is at y = 0). 
Since the Green function decays as l/x 2 , we can only 
have n terms off the block diagonal. Furthermore, to get 
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a connected correlation function, we must have exactly 
one term off the block diagonal in every block row and in 
every block column. This allows us to generalize Eqs. (31- 
32) for open boundary n-point functions; they generalize 
in the obvious manner, with (n — 1)! trace terms for the 
n-point function, corresponding to the (n — 1)! ways that 
we can loop through the n positions. 

Eq. (33) shows that each off-diagonal block, G uv , fac- 
torizes into the product of a column vector and row vec- 
tor: 



G 



1 



Ufl 



rh„,hf, 



(34) 



1 of the sites 



where h„ , is a column vector of heights y - 
around Xf in L oPiUf (xf) — i.e. the p th entry of h Uf is the 
value of y + 1 for the p th site of L oPiU (xf). Substituting 
this in the generalization of Eqs. (31-32), and using the 
cyclicity of the trace, each of the (n — 1)! matrix traces 
becomes a product of n 1 x 1 matrices. Furthermore, 
the [n — 1)! traces differ from each other only in the 
l/(xf — x g ) 2 terms chosen. The leading-order, connected 
part of the correlation function of n local arrow diagrams 
is then found to be 



(Lop,u 1 (xi)L op ,u 2 (x2) ■ ■ ■ L op . Un (x n )} 



= ( n ) dct m ( 35 ) 



M is defined as the n x n matrix 



M 



r o iff= g 

f9 I y(Xf Xg) 2 if f ± g 



(36) 



and the k u are simply numbers: 



Pa, op is the probability for a site along an open boundary 
to have height a, and the K a are normalization factors. 
We then, finally, have 



(0oi(a;i)</»o2 (x 2 ) 



,(a; n ))=det(M) (40) 



For n — 2, this reproduces the open boundary one- and 
two-point functions found in [16]. 

det(M) is the same as the n-point function of —28889, 
so up to rescaling factors (— 2if a 's), all four height vari- 
ables are represented by 8989 along open boundaries. 
This is rather surprising, given that we have seen that 
the height variables are represented by different operators 
along closed boundaries (Eqs. (20-22)). InCFT's, bound- 
ary operators are derived from OPE's of bulk operators — 
so the fact that the height operators are different along 
closed boundaries proves that they must be different in 
the bulk, but apparently these different bulk operators 
become identical along open boundaries. 

We nowhere used the fact that these were the local 
arrow diagrams associated with the height variables. So, 
in fact, we have shown that all local arrow diagrams along 
open boundaries are represented by 8989. 

We have also found the correlation function of n unit 
height variables along closed boundaries. This requires 
local arrow constraints at 3n vertices of the ASM, and 
thus the calculation of a 3n-dimcnsional matrix determi- 
nant. The matrix is divided into 3 by 3 block submatri- 
ces, such that the diagonal blocks are all identical, and 
the off-diagonal blocks all have the same form. A rotation 
makes the matrix diagonal in 2 out of every 3 rows (and 
columns). The universal part of the correlation function 
is thus found to be 



37T-8 



dctM 



(41) 



k u = - dct (I + B„G„„) (hi — i— 



-B„h, 



(37) 



Inserting this into Eq. (30) gives all open boundary 
n-point correlations. To express our results in a simpler 
manner, we define 



4>a{x) 



Sh x ,a - Pa 



op 



.where a =1,... 4 (38) 



We have defined the following constants: 



Pl,op = 


9 
2 


42 + 


320 


512 
9tt 3 


K x 


P2,op = 


33 
4 


, 66 

IT 


_ 160 
tt' 2 


, 1024 


K 2 


P3,op = 


15 

4 


22 , 

TV 


160 
Ztt 1 


512 

" 97?^ 


K 3 


P4,op = 


1 - 


2 

7T 






K 4 



_ 3 , 80 
tt 3-tt 2 9tt 



512 



TT 3-7T 2 



1024 
9tt 3 



7_ , 40 _ 512 

TV ~TT^ 9TT 3 



(39) 



This confirms the field identification in Eq. (20). 



XI. n-POINT CORRELATIONS ALONG OPEN 
BOUNDARIES, WITH DISSIPATIVE DEFECTS 

Along open boundaries, the defect-free Green function, 
G = G OP:0 docs not diverge as L — > oo, so for a single 
dissipative defect we can modify the Green function as 
in Eq. (25). Using this new Green function, the open 
height probabilities at (xi,0), for a defect of strength k 
at d = (0, y) are 



f op (a) = -K a 



k(y+l) 2 



7r(l + kG oPt0 (d,d)) 



a = 1,2,3,4 
(42) 

The same K a factors that we saw in the height-height 
correlations appear in height-defect correlations. 
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We define an operator (f>5-k(d), corresponding to the 
addition of a defect of strength k at d — (x, y), and then 
multiplication of all correlation functions by a normal- 
ization factor 

tt(1 + kG opfi (d,d)) 

k( y + iy (4dJ 

Then Eq. (42) becomes 

(<t>a{xi)<f>6;k(X2)) = —. ~ 7J, a^h (44) 

\X\ %2) 

05 ; fc acts just like any of the four height variables in two- 
point correlations (Eq. (44) is Eq. (40) with n = 2). In 
fact, we find that </>s ; fc acts like <fii, <j>2, 4>z, and </>4 in 
all higher-order correlation functions, containing multiple 
height variables and multiple dissipativc defects. 

Suppose we are calculating a correlation function with 
n height variables, and m dissipative defects. The dissi- 
pative defects are at d w = (x w ,y w ), and have strength 
k w , 1 < w < m. As with the height locations, the x w 
coordinates of the defects all scale with the same factor 
x, where x — > 00. The change in the toppling matrix , 
5 A = A - A , is 

j A ,, = f k w if i =.f = d w , l<w<m ^ 
\ otherwise 

The Green function is modified from its defect-free value, 
G p,o, to 



(We already saw this factor for a single dissipative defect, 
in Eq. (25).) Furthermore, inspection of Eq. (33) shows 

that the visit to the defect at d w from another site will 
result in a factor of k w (y w + l) 2 / 71 "- With Eq. (47), this 
motivates the normalization factor in Eq. (43). 

Eq. (43) normalizes the correlation function of n height 
variables and m defects. To see that the form of the cor- 
relation function is still det(M), expand the determinant 
out into cycles. The connected part of the determinant 
in Eq. (40) is a sum of closed cycles of length n, where 
each cycle visits each of the positions (xf , 0) exactly once, 
and picks up a factor of l/(x/ — x g ) 2 when it travels from 
xj to x g . After normalizing, the defects have exactly the 
same effect as the height variables — each trip to (or from) 
a defect results in a 1/x 2 term from the Green function 
to (or from) the defect (Eq. (33)). 

So, in the end, the correlation function of n height 
variables on the boundary and m defect sites near or 
on the boundary, is given by the (m + n) dimensional 
matrix determinant, det(M) (with appropriate normal- 
ization factors). This shows that dissipative defect sites 
along or near open boundaries are, like the height vari- 
ables, represented by 8989. 

Note that a dissipative defect has a much larger ef- 
fect along a closed boundary than along an open one. 
A defect is represented by a dimension zero operator 
along a closed boundary, but by a dimension two op- 
erator along an open boundary. This makes sense; along 
open boundaries, grains are already dissipated by top- 
plings, so adding a little more dissipation has only minor 
effects, compared to dissipation on a closed boundary. 



G 1 - 1 - Gop '° 

A A Q + SA I+(<5A)G op , 

OO 

= ^G op ,o(-(<5A)G op ,o) p (46) 

G(i,j) can be represented as a trip from i to j, where 
along the trip, the traveller can visit any of the defect 
sites as often as he or she wishes, each time picking up a 
factor of — (5A)G op . . 

We have already seen that the defect-free correlation 
function of n height variables has a leading term of 
0(l/x 2n ). If we instead use the Green function with de- 
fects, each trip to a defect introduces a factor of 1/x 2 (sec 
Eq. (33)). In a connected function, we should visit each 
defect at least once; in the leading term, each defect will 
be visited from a distant site exactly once, and the corre- 
lation function will have a leading term of 0(l/x 2 ^ n+m ^). 

After visiting d w , we may travel repeatedly from d w 
to d w without picking up extra factors of 1/x 2 . This 
produces a contribution to Eq. (46) of 

00 - 1 

^(-^0(4,4))^ = ——^ (47) 
p= o l + k w Uo(d w , a w ) 
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APPENDIX A: ANOMALOUS GRAPHS IN 
BOUNDARY TWO-POINT 
CORRELATIONS — PART I 

In this section we discuss what we call anomalous 
graphs of the first kind, which arise when converting from 
ASM height probabilities to spanning tree probabilities. 
As stated in section V, it would be natural to expect, 
based on analogy with the one-point height probabili- 
ties, for the two-point height probabilities to be given 
by Eq. (5). However, this turns out to not be the case. 
Let us carefully consider how height correlations can be 
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FIG. 7: State not in any S^, and in multiple Ski- 




FIG. 8: Venn diagram in the space of states where (/i?, hj) — 
(2, 2) is allowed, but (ft.?, hj) — (1, 1) is not. 

turned into spanning tree probabilities. We focus on the 
closed two-point correlations; other cases are similar. 

For correlations between i and j, Ivashkevich divided 
the states of the ASM into sets Ski , consisting of states 
allowed when ft? > fc and hj > I, but forbidden oth- 
erwise [16]. However, not all ASM states fall into one 
of these sets. There are states that are allowed when 
(ft?, ft?) = (1,2), and when (ft?, ft?) = (2,1), but forbid- 
den when (ft?, ft?) = (1,1); these states do not belong to 
any set Ski- One such state is shown in figure 7. 

We find it convenient to define Ski , consisting of ASM 
height configurations on the sandpile, excluding i and j, 
that are allowed when we add (ft?, ft?) = (k,l), but both 
forbidden when we add (ft?, ft?) = (k — 1,1), and also 
forbidden when we add (ft?, ft?) = (k,l — 1). Note that 
when counting states, the fact that we do not specify the 
heights of i and j introduces a multiplicative factor; for 
example, | \ — \ S12 | /(3 * 2) for i and j both on a 
closed boundary. Now, every state must be in at least 
one of the sets Ski, but some states are in several Ski's. 
For example, the state in figure 7 is in both S\2 and S21 ■ 

No anomalous graphs arise for two-point correlations 
involving at least one unit height variable, since in those 



cases the relevant S'-sets do not intersect. The number 
of states where i has height one and j has height ft is 
S p =i I Sip I- Furthermore, the representation of S\h 
is exactly what we would expect; it corresponds to the 
number of spanning trees where NNP? = 0, and NNP? = 
h — 1 . So no anomalous graphs arise when the two-point 
correlation has at least one unit height variable. 

Things get more complicated when both heights are 
higher heights. We discuss in detail the 2-2 (height two- 
height two) correlation along a closed boundary; the anal- 
ysis for the other two-point correlations is similar. 

If i and j both have height two, we must be in at least 
one Ski, for k < 2, I < 2. S12 and S21 intersect, so 
number of 2-2 states is 



I Sn I + I S12 1 + I S21 I + (l S 22 1 - I S 12 n S21 |) 

(Al) 

The first three terms all have the "natural" span- 
ning tree representation. The difficulty is in evaluat- 
ing I S22 I — I S12 n S21 |. 5*22 counts states where 
(ft?, ft?) = (2,2) is allowed, but neither (ft?, ft?) = (1,2) 
nor (ft?, ft?) = (2, 1) are allowed. In figure 8, wc have 
started with a large rectangle, representing the the set of 
states where (ft?, ft?) = (2, 2) is allowed, and (ft?, ft?) = 
(1,1) is forbidden — we call this set X. In the rect- 
angle are two subsets, corresponding to regions where 
(ft?, ft?) = (1,2) is allowed, and where (ft?, ft?) = (2, 1) is 
allowed — we call these two sets X\ and X 2 . In this Venn 
diagram, 522 is the diagonally shaded region outside the 
circles, and S12 (~l S 2 i is the cross-hatched intersection of 
the two circles. Looking at the Venn diagram, we see 
that to find | S22 | — | §12 H S21 |, we start with the set 
X, and then subtract off the states in X\ and X2 inde- 
pendently. By independently we mean that states in the 
intersection of X\ and X 2 get subtracted off twice. So 

I §22 I - I S12 nS 2 i\=\X\-\X 1 \-\X 2 \ (A2) 

For states in X, setting (ft?, ft?) = (1,1) produces an 
MFSC. (Note that we have defined the MFSC as the 
maximal FSC produced when the heights at both i and 
j are simultaneously set to one; if only one height was 
set to one, then the largest FSC might be smaller, or 
there might be no FSC at all.) The set X can be parti- 
tioned into the following disjoint subsets, depending on 
the shape of the MFSC: 



Xa ■ The MFSC consists of disjoint subsets around i and j 

Xb ■ The MFSC is connected, with > 2 neighbors of i and exactly 1 neighbor of j 
Xc ■ The MFSC is connected, with > 2 neighbors of j and exactly 1 neighbor of i 
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Xd ■ The MFSC is connected, with exactly 1 neighbor of i and exactly 1 neighbor of j 



X\ and X 2 can be partitioned into analogously defined 
subsets, Xia, Xib, etc. . . So, for example, X\u is the 
subset of Xjj such that (/ij, hj) = (1, 2) is allowed. (Note 
that Xic = and X 2 b = 0.) We want 



E 



Xk 



Xu 



X 2k |) 



(A3) 



k£{A,B,C,D} 



Our "naive" guess would be that this would equal the set 
of states where NNPt- = NNPj = 1 . We carefully count 
the states, comparing with this guess. 

Xa, after subtracting off the states from Xia and X 2 a, 
is equal to one-fourth the number of spanning trees for 
which i and j each have one NNP, and {i}UTree--does not 
border or intersect {j} U Treej-. (The one-fourth comes 

from the fact that the spanning tree arrows from i and 
j can point out from the MFSC in any direction.) Treej- 
refers to the set of sites that are predecessors of i. The 
condition that {«} U Treer and {j} U Treej- cannot bor- 
der each other comes from the condition that the MFSC 
consist of disjoint subsets around i and j. 

We now consider X B . The MFSC generated when 
(%, hj) = (1, 1) must still be an MFSC when (%, hj) = 
(2,1). So (%,/ij) = (2,2) (2,1) produces an MFSC 

that includes exactly one neighbor of j, the site i, and 
at least two of i's neighbors. Just as in section III, this 
is equivalent to a modified ASM, where bonds along the 
border of the MFSC are removed (except for one bond 
of j). In Eq. (A3), X 2 b = 0, but we do need to subtract 
off states in Xib- To do this we only count the states of 
Xb such that (/ij, hj) = (1,2) is forbidden, which implies 
that /ij = 2 — ► 1 should produce a new, smaller FSC, 
completely contained within the larger MFSC. We then 
see that \ Xb \ — \ X\ B \ — \X 2 b | corresponds to one- 
fourth of all spanning trees where NNPt- = NNPj* = 1, 

and either i e Treej- or i borders Treej-, with one ex- 
ception. The exception occurs because the MFSC, by 
construction of Xb, can only have one neighbor of j. 
So | Xb I — | X\b | — | X 2 b | will not count cases 

where NNPt- = NNPt- = 1 , i borders Trccr, and j bor- 
i j j 

ders Tree^ — this case is shown in figure 3. We label this 
set of graphs as c. Since c has NNPt- = NNPj- = 1, it 
would be natural to expect it to appear in the spanning 
trees contributing to the height two-height two correla- 
tion function. However, since no MFSC's of X have two 
neighbors of i and two neighbors of j, graph c does not 
appear in Xc or Xd either. 

The analysis for Xq is identical to that for Xb, and 
Xc | — | Xic | — | X 2 q | counts one-fourth the span- 
ning trees where NNPt- = NNPt- = 1 , and either j £ Treey 

or j borders Treej-, except that, again, the spanning trees 
of c are excluded. 




1 J 

FIG. 9: Correlations where fi\ / fii 



In Xb, the MFSC has only one neighbor of i and one 
neighbor of j. The one-to-one mapping between Xb and 
spanning tree states posessess some subtleties, but the 
end result is what one would expect: The number of 
states in | Xd \ — \ Xm | — | X 2 d \ is one-fourth the 
number of spanning trees where NNPt- = NNPt- = 1 , and 

Treer and Treej- border each other, but i £ Treer, and 

j £ Tree^. 

In the end, we see that 4 ^| S 22 \ — \ S\ 2 (~l 5*21 (j is 
equal to the number of spanning trees where NNPt- = 
NNPt- = 1 , except for the set c, which contains all span- 
ning trees where both i borders Treej-, and j borders 
Treej-. c consists of the anomalous graphs of the first 
kind. 

An similar analysis for other closed correlation func- 
tions shows that the spanning trees in c contribute to the 
height two-height three correlation, and get subtracted 
from the height three-height correlation (relative to a 
"naive" approach). These changes are necessary for the 
height probabilities to all sum to one, so this provides 
a check on our mapping between ASM states and span- 
ning tree states. In the next appendix, we consider yet 
another complication that arises in the calculation of the 
two-point functions. 



APPENDIX B: ANOMALOUS GRAPHS IN 
BOUNDARY TWO-POINT 
CORRELATIONS— PART II 

We saw in the previous appendix that Eq. (5) does 
not quite hold, but is only off by the anomalous graph 
c. So except for this complication, the two-point height 
correlations can be turned into linear combinations of 
probabilities for spanning trees with nonlocal conditions 
on NNPt- and NNPj. As in the previous appendix, we dis- 
cuss only the closed boundary two-point functions; other 
cases are similar. 

As discussed in section IV, we can write the probability 
Af h to have NNPt- = h - 1 linear combination of 

nonlocal arrow diagrams, which we can then rewrite as 
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b term 



FIG. 10: Correlations where <f>\ / <f>2 



follows (see figures 1 and 2): 

Af 2 = 2(/3 1 +/3 2 + 7 ) 

= 2(-L Cil + 3L Ci2 + L Ci3 ) + 

4(/3i - 02) + 2(0i - <h) (Bl) 
A/" 3 = 2(0i + 2 + e) + * 

= 2(2L c ^i — 2L C ^2 + L c ^) + L Ct 4 — 

4(/3i - 2 ) ~ 2(0i - 2 ) (B2) 

These relationships hold regardless of the correlation 
functions that i are in. We define operators correspond- 
ing to the local parts of these terms: 

Ln2 = 2(— L C; i + 3L C ,2 + ^0,3) (B3) 
L N3 = 2(2L cA - 2L C , 2 + L c ,5) + £ c ,4 (B4) 

In one-point functions, /3i = /?2, and 0i =02, so A/2 
and A/3 simply become Ljv 2 and £jv3, which are lo- 
cal, and whose expectation values can be found with 
the Majumdar-Dhar method, as discussed in section IV. 
(These are the same relations found in [16].) 



It would be simplest if in correlations of A/2 and A/3, 
we could replace A/2 and A3 with L^2 and Ljq 3 , since 
local correlation functions are easily calclated. As in the 
previous appendix, we call this approach "naive," — this 
naive approach does not quite work, and we call the de- 
viations of the correct answers from the naive approach 
"anomalous graphs of the second kind." 

We no longer have 0\ = 02 and 0i = 02 in correlations 
of A/2 and A3, because in correlations between distant 
sites i and j, switching arrows at i can affect predeces- 
sor relationships at j. To analyze the (0\ — 02) terms 
in Eqs. (B1-B2), consider the configuration in figure 9, 
where we have not specified the direction of the arrow 
from i. If the arrow from i points to i\, then ji is not 
a predecessor of J, so the configuration at j is 7. And if 
the arrow from i points to 12, then j\ is a predecessor of 
j, so the configuration at j is 2 . So switching from (5\ to 
02 at i can affect whether the configuration at j is 02 or 
7. However, this inequivalence between fj\ and 02 turns 
out to have no effect on any correlation functions, to any 
order, since 2 and 7 always appear in the combination 
02 + 7, in Ln2 and Ln 3 , and 02 + 7 has no correlations 
with 0i —02- 



Unfortunately, things become more complicated with 
the 0i and 02 terms. If the long path of 0i avoids arrow 
restrictions at other sites, the long path can be reversed, 
and 0i will be equivalent to 2 . However, if the long 
path goes through arrow restrictions at other sites of the 
correlation function, then 0! will not be equivalent to 02- 
Figure 10 compares diagrams that arise in (0i, 0i) corre- 
lations, with diagrams that arise in (02, 02) correlations. 
Three of the four diagrams shown (labeled with "a" 's) 
are equivalent, but the fourth one (labeled with "b") is 
not. The resultant anomalous graphs of the second kind 
were shown and discussed earlier, in section V (see fig- 
ure 4). When all the correlations of the (0i — 02) 's in 
Eqs. (B1-B2) are considered, we find the following for 



closed boundary correlations: 



<A/2(r)A/2(0)> = (L N 2(r)L N2 (0))+4(a-b) (B5) 
(AA 2 (r)AA 3 (0)) - (L N2 (r)L N3 (0)} -A(a-b) (B6) 
(M 3 (r)M 3 (0)} - (L N3 (r)L N3 (0))+A(a-b) (B7) 

The correlations involving the height one variable are un- 
affected by these complications. 

To summarize, naively transforming from nonlocal ar- 
row diagrams to local arrow diagrams independently at 
every site of a two-point correlation results in anomalous 
graphs a and b. 
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APPENDIX C: GREEN FUNCTIONS 



The inverse of the bulk toppling matrix A is the lat- 
tice Green function, which has long been known [24]. It 
is given by 



o[x,y > J Q 2tt i 2?r 4 - 2cospi - 2cosp 2 

(CI) 

This integral is divergent, producing terms of order InL, 
where L is the system size, but these divergences are usu- 
ally unimportant, since we are typically concerned with 
differences in Green functions. For large x, y this has the 
expansion [24] 



1 



G (x, y) = -— \n(x 2 + y 2 ) - J- - =p + 



In 8 



4ir 



2tt 4ir 



(C2) 



where 7 = 0.57721 ... is the Euler-Mascheroni constant. 
For sites (xi,yi) and (2:2,2/2) near an open boundary, 
where x is the coordinate along the boundary, and y is 
the distance from the boundary (located at y = 0), the 
Green function is [25] 



G op ,o(x 1 ,y 1 ;x 2 ,y2) = G (xi - x 2 ,yi - 2/2) - 

G (x 1 -x 2 ,y 1 + y 2 + 2) (C3) 

Along closed boundaries, it is [25] 

G c ,o(zi,2/i; £2,2/2) = Go (#i - #2,2/1 - 2/2) + 

G (xi- x 2 , yi +y 2 + l) (C4) 

The minus sign between the Green functions in Eq. (C3) 
cancels out divergences in the Green function. The ex- 
pansion of the Green function for points along the bound- 
ary has already been calculated [25], and can be extended 
to points near the boundary, but far from each other 
(2/1 = 0(1), y 2 = 0(1), I X! - x 2 |-> 00), by the recur- 
sion relationships, GqA = I. We find 



G op {xi, yi; x 2 ,y 2 ) = 



ir(xi - x 2 ) 2 



+ ... 



(C5) 



and 



J 



G c (x 1 ,y 1 ,x 2 ,y 2 ) = -- In | x Y - x 2 \ -(- + ^^) - (3j/i(j/i + 1) + 3y 2 (y2 + 1) + 1)7-7 — " ^2 

7T 7T 27T 07T(Xl — X 2 ) 
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+ 2/1 (yi + l)(2/i + 2/1 - 1) + 2/2(2/2 + l)(2/ 2 + 2/2 - 1) + 62/1(2/1 + 1)2/2(2/2 + 1) - 77; 



1 



60 y 4tt(xi - x 2 ) 



+ ...(C6) 



APPENDIX D: MORE CLOSED BOUNDARY 
CORRELATION FUNCTIONS 

Here we list the three-point correlation functions along 
closed boundaries that were not stated in section VI. As 
a check, the correlation functions in this appendix were 
found by the methods already described in sections V 
and VI. However, they can all also be determined from 
those already listed in section VI, from the requirement 
that the three height probabilities must sum at all sites, 
and by symmetry. They are listed here only for reference, 
and because they provide checks on our calculations. We 
find 



(3tt-. 



+ 



/c(l,3,3) = 



7T 5 (X1 - X 2 ) 3 (x 2 - X 3 ) 3 

(3tt - 8) (24- 7tt) 
27r 5 (xi - a; 2 ) 3 (a;i - x 3 ) 3 
(37T-8)(4-7r)(8 + 7r) 
tt 6 (xi - a; 2 ) 2 (a;i - x 3 ) 2 (x 2 - x 3 ) 
(37r-8)(8-7r) | ^ 
2tt 5 (xi — x 2 ) 3 (xi — x 3 ) 3 



+ 



(D2) 



(D3) 



We can now check that / c (l, 1, 1) + / c (l, 1,2) + 
f c (l,l,3) = 0, as it must. Interchanging x 2 and X3 in 
/ c (l,l,2) gives 



/c(l,l,3) = - 



/c(l,2,3) = 



2(4- 7r)(37r 



+ 



7r 6 (a;i - x 2 ) 2 (xi - x 3 ) 2 (x 2 - x 3 ) 2 
{'in - 8) 2 



+ ... 



7T 5 (X1 - X 3 ) 3 (X 2 - X 3 ) 3 

4(7r-3)(7r + 8)(37r-8) 
7r 6 (a;i - x 2 ) 2 (xi - x 3 ) 2 (x 2 - x 3 ) 2 



(Dl) 



/c(l,2,l) = - 



8(7T-3)(37T- 



7T 6 (X1 - X 2 ) 2 (X 1 - X 3 ) 2 {X 2 - X 3 ) 2 
(37T-8) 2 



7T 5 (X1 — X2) 3 (x 3 - X 2 ) 3 



+ ... 



(D4) 



We can then check that / c (l,l,l) + / c (l,2,l) + 
/ c (l, 3, 1) = 0. Three-point correlation functions with no 
unit height variables cannot be found with the methods in 
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this paper, as already discussed in section VI. However, if 
we use the conjecture proposed in section V, of dropping 
the anomalous graphs (as we did to obtain Eq. (23)), we 
now obtain 



/c(2,2,3) 



/c(2,3,3) 



(24 - 5tt)(-192 + 112tt - 13tt 2 ) 

47T 6 (X1 - X 2 ) 2 (X1 - X 3 ) 2 (X 2 - X 3 ) 2 

(37r-8)(77r-24) + 

27T 5 (xi — X 3 ) 3 (x 2 — X 3 ) 3 

(8 + tt)(-192 + 112tt- 13tt 2 ) 
4it 6 (xi - x 2 ) 2 {xi - x 3 ) 2 (x 2 - x 3 ) 2 
(37r-8)(8-7r) 



+ . . . 



/ c (3,3,3) = - 



2ir 5 (xi — x 2 ) 3 (xi — X3) 3 

(8 + 7r)(64 - 327r + 7T 2 ) 
47r 6 (xi - x 2 ) 2 (x 1 - x 3 ) 2 (x 2 - x 3 ) 2 



(D5) 



(D6) 



+ 



(D7) 



As with the other correlation functions in this section, 
these agree with the requirements that the three height 
probabilities must sum to one at all sites, and with the 
field identifications in Eqs. (20-22). 



APPENDIX E: THE c = -2 CONFORMAL FIELD 
THEORY 

The central charge -2 conformal field theory is per- 
haps the simplest known logarithmic conformal field the- 
ory. While the theory has a simple underlying Gaussian 
structure, it still possesses a number of subtle features. 
We use the formulation of the c = — 2 CFT where the 
action is given by 



should still be expected to apply. In particular, just as 
for non-logarithmic boundary CFT's, as bulk operators 
are moved near a boundary, their antiholomorphic pieces 
should behave like holomorphic pieces at mirror locations 
across the boundary [19, 30]. 
APPENDIX F: A PROOF THAT THE HEIGHT 
VARIABLES HAVE DIFFERENT BULK FIELD 
IDENTIFICATIONS 

The correlation functions in Eqs. (13-15), Eqs. (17-19), 
and Eqs. (27-29) show conclusively that the three height 
variables are represented by different operators along 
closed boundaries. As already discussed, since bound- 
ary operators are derived from OPE's of bulk operators, 
this proves that the height variables must be represented 
by different operators in the bulk as well [19]. However, 
it is worth noting that this conclusion can be reached 
with a simple argument, based on general principles of 
conformal field theory, without doing any detailed calcu- 
lations. 

Suppose that all four height variables were represented 
(up to multiplicative factors) by the same field operator. 
The unit height variable is known, from its two-point 
correlation, to have dimension two [14], so, by our as- 
sumption, all four height variables would have scaling di- 
mension two. The height probabilities get modified from 
their bulk values, ps,h {h = 1,2,3,4) near a boundary 
(closed or open). Then one-point functions of operators 
of dimension d will decay as l/y d , where y is the dis- 
tance from the boundary, and d is the operator dimen- 
sion [28, 30]. 



Ph{y) =PB,h 



Ch 
„2 



(Fl) 



s = \ j deoe (ei) 

d and d refer to the holomorphic and antiholomorphic 
derivatives — that is, the derivatives with respect to z = 
x + iy and z = x — iy. and 9 are anticommuting Grass- 
manian variables. The action has zero modes, which 
make the partition function zero. If we normalize the 
action by not integrating over the zero modes, we get 
Wick contraction rules, with each contraction between 
0(z\) and 0(z 2 ) giving a factor of — log(zi — z 2 ). 

While boundary conformal field theories are generally 
well understood [19, 28], boundary logarithmic confor- 
mal field theories possess a number of subtleties that 
are not well understood. Results on boundary LCFT 
are still to some degree contradictory [29-34]. How- 
ever, basic results from non-logarithmic boundary CFT 



for some constants c^. If the fields arc normalized (to 
have zero expectation value and coeficient -1 in two-point 
correlations), then general CFT principles predict that 
the coefficients of the 1/y 2 terms should be universal 
numbers, depending only on the field and the boundary 
condition [28]; in particular, they should be independent 
of h. So upon normalizing the height variables, the dif- 
ferent Ch should all become c, a number independent of 
h. Since we are assuming that all four height variables 
are represented by the same field, the 1-1, 2-2, 3-3, and 
4-4 correlations should all have the same sign (negative), 
so this normalization should not change the signs of the 
coefficients, and all the c^'s should have the same sign as 
c. However, we need c h = 0> f° r the f° ur height 

probabilites to sum to one, so the Ch cannot all have the 
same sign. By contradiction, the four height variables 
must be represented by different fields in the bulk. 
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